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Abstract 

We present C programming language versions of earlier published Fortran programs (Muruganandam 
and Adhikari (2009) [1]) for calculating both stationary and non-stationary solutions of the time-dependent 
Gross-Pitaevskii (GP) equation. The GP equation describes the properties of dilute Bose-Einstein con- 
densates at ultra-cold temperatures. C versions of programs use the same algorithms as the Fortran ones, 
involving real- and imaginary-time propagation based on a split-step Crank-Nicolson method. In a one-space- 
variable form of the GP equation, we consider the one-dimensional, two-dimensional, circularly-symmetric, 
and the three-dimensional spherically-symmetric harmonic-oscillator traps. In the two-space-variable form, 
we consider the GP equation in two-dimensional anisotropic and three-dimensional axially-symmetric traps. 
The fully-anisotropic three-dimensional GP equation is also considered. In addition to these twelve programs, 
for six algorithms that involve two and three space variables, we have also developed threaded (OpenMP 
parallelized) programs, which allow numerical simulations to use all available CPU cores on a computer. All 
18 programs are optimized and accompanied by makefiles for several popular C compilers. We present typ- 
ical results for scalability of threaded codes and demonstrate almost linear speedup obtained with the new 
programs, allowing a decrease in execution times by an order of magnitude on modern multi-core computers. 
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New version program summary 

Program title: GP-SCL package, consisting of: (i) imagtimeld, (ii) imagtime2d, (iii) imagtime2d-th, (iv) imagtimecir, 
(v) imagtime3d, (vi) imagtime3d-th, (vii) imagtimeaxial, (viii) imagtimeaxial-th, (ix) imagtimesph, (x) realtimeld, 
(xi) realtime2d, (xii) realtime2d-th, (xiii) realtimecir, (xiv) realtime3d, (xv) realtime3d-th, (xvi) realtimeaxial, (xvii) 
realtimeaxial-th, (xviii) realtimesph 

Journal Reference: Comput. Phys. Commun. 183 (2012) 2021. 
Catalogue identifier: AEDU_v2_0 

Program Summary Lfi?£/ [http://c pc.cs.qub.ac.uk/summaries/AEDU_v2_0.htmll 
Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland 
Licensing provisions: Standard CPC licence, [http://cpc.cs.qub.ac.uk/lice nce/licence.html 
No. of lines in distributed program, including test data, etc. : 180 583 
No. of bytes in distributed program, including test data, etc.: 1 188 688 
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Distribution format: tar.gz 

Programming language: C and C/OpenMP 

Computer: Any modern computer with C language compiler installed 
Operating system: Linux, Unix, Mac OS, Windows 

RAM: Memory used with the supplied input files: 2-4 MByte (i, iv, ix, x, xiii, xvi, xvii, xviii), 8 MByte (xi, xii), 32 
MByte (vii, viii), 80 MByte (ii, iii), 700 MByte (xiv, xv), 1.2 GByte (v, vi) 

Number of processors used: For threaded (OpenMP parallelized) programs, all available CPU cores on the computer. 
Classification: 2.9, 4.3, 4.12 

Catalogue identifier of previous version: AEDU_vl_0 

Journal reference of previous version: Comput. Phys. Commun. 180 (2009) 1888. 
Does the new version supersede the previous version?: No 

Nature of problem: These programs are designed to solve the time-dependent Gross-Pitaevskii (GP) nonlinear par- 
tial differential equation in one-, two- or three-space dimensions with a harmonic, circularly-symmetric, spherically- 
symmetric, axially-symmetric or fully anisotropic trap. The GP equation describes the properties of a dilute trapped 
Bose-Einstein condensate. 

Solution method: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretizing 
in space and time. The discretized equation is then solved by propagation, in either imaginary or real time, over 
small time steps. The method yields solutions of stationary and/or non-stationary problems. 

Reasons for the new version: Previous Fortran programs [1] are used within the ultra-cold atoms [2-11] and nonlinear 
optics [12,13] communities, as well as in various other fields [14-16]. This new version represents translation of all 
programs to the C programming language, which will make it accessible to the wider parts of the corresponding 
communities. It is well known that numerical simulations of the GP equation in highly experimentally relevant 
geometries with two or three space variables are computationally very demanding, which presents an obstacle in 
detailed numerical studies of such systems. For this reason, we have developed threaded (OpenMP parallelized) ver- 
sions of programs imagtime2d, imagtime3d, imagtimeaxial, realtime2d, realtime3d, realtimeaxial, which are named 
imagtime2d-th, imagtime3d-th, imagtimeaxial-th, realtime2d-th, realtime3d-th, realtimeaxial-th, respectively. Fig- 
ure 1 shows the scalability results obtained for OpenMP versions of programs realtime2d and realtime3d. As we 
can see, the speedup is almost linear, and on a computer with the total of 8 CPU cores we observe in Fig. 1(a) a 
maximal speedup of around 7, or roughly 90% of the ideal speedup, while on a computer with 12 CPU cores we 
find in Fig. 1(b) that the maximal speedup is around 9.6, or 80% of the ideal speedup. Such a speedup represents 
significant improvement in the performance. 




number of CPU cores number of CPU cores 

Figure 1: (Colour online) Speedup in the execution time of realtime2d-th and realtime3d-th threaded (OpenMP parallelized) 
programs as a function of the number of CPU cores used. The results are obtained: (a) on an 8-core machine with 2 x quad-core 
Intel Nehalem E5540 CPU at 2.53 GHz, using the ice compiler, (b) on a 12-core machine with 2 x six-core Intel Nehalem X5650 
CPU at 2.66 GHz, using the pgee compiler. The spatial grid sizes used are 2000x2000 (realtime2d-th) and 1000x1000x300 
(realtime3d-th) . 



Summary of revisions: All Fortran programs from the previous version [1] are translated to C and named in the 
same way. The structure of all programs is identical. We have introduced the use of comprehensive input files, 
where all parameters are explained in detail and can be set by a user. We have also included makefiles with tested 
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and verified settings for GNU's gcc compiler, Intel's ice compiler, IBM's xlc compiler, PGFs pgee compiler, and 
Oracle's suncc (former Sun's) compiler. In addition to this, 6 new threaded (OpenMP parallelized) programs are 
supplied (imagtime2d-th, imagtime3d-th, imagtimeaxial-th, realtime2d-th, realtime3d-th, realtimeaxial-th) for algo- 
rithms involving two or three space variables. They are written by OpenMP-parallelizing the most computationally 
demanding loops in functions performing time evolution (calcnu, calclux, calcluy, calcluz), normalization (calcnorm), 
and calculation of physical quantities (calcmuen, calcrms). Since some of the dynamically allocated array variables 
are used within such loops, they had to be made private for each thread. This was done by allocating matrices 
instead of arrays, with the first index in all such matrices corresponding to a thread number. 

Additional comments: This package consists of 18 programs, see Program title above, out of which 12 programs (i, 
ii, iv, v, vii, ix, x, xi, xiii, xiv, xvi, xviii) are serial, while 6 programs (iii, vi, viii, xii, xv, xvii) are threaded (OpenMP 
parallelized). For the particular purpose of each program, please see descriptions below. 

Running time: All running times given in descriptions below refer to programs compiled with gcc on quad-core Intel 
Xeon X5460 at 3.16 GHz (CPU1), and programs compiled with ice on quad-core Intel Nehalem E5540 at 2.53 GHz 
(CPU2). With the supplied input files, running times on CPU1 are: 5 minutes (i, iv, ix, xii, xiii, xvii, xviii), 10 
minutes (viii, xvi), 15 minutes (iii, x, xi), 30 minutes (ii, vi, vii), 2 hours (v), 4 hours (xv), 15 hours (xiv). On CPU2, 
running times are: 5 minutes (i, iii, iv, viii, ix, xii, xiii, xvi, xvii, xviii), 10 minutes (vi, x, xi), 20 minutes (ii, vii), 1 
hour (v), 2 hours (xv), 12 hours (xiv). 

New version program summary (i) 

Program title: imagtimeld 

Title of electronic files: imagtimeld. c, imagtimeld. h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 4 MByte 

Programming language used: C 

Typical running time: 2 minutes (CPU1), 1 minute (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in one space dimension with a harmonic trap. The GP equation describes the properties of a dilute trapped 
Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in imaginary time over small time 
steps. The method yields solutions of stationary problems. 

New version program summary (ii) 

Program title: imagtime2d 

Title of electronic files: imagtime2d.c, imagtime2d.h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 80 MByte 

Programming language used: C 

Typical running time: 30 minutes (CPU1), 20 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in two space dimensions with an anisotropic trap. The GP equation describes the properties of a dilute 
trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in imaginary time over small time 
steps. The method yields solutions of stationary problems. 

New version program summary (iii) 

Program title: imagtime2d-th 

Title of electronic files: imagtime2d-th.c, imagtime2d-th.h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 80 MByte 

Programming language used: C/OpcnMP 

Typical running time: 15 minutes (CPU1), 5 minutes (CPU2) 
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Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in two space dimensions with an anisotropic trap. The GP equation describes the properties of a dilute 
trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in imaginary time over small time 
steps. The method yields solutions of stationary problems. 

New version program summary (iv) 

Program, title: imagtimecir 

Title of electronic files: imagtimecir. c, imagtimecir. h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 2 MByte 

Programming language used: C 

Typical running time: 2 minutes (CPU1), 1 minute (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in two space dimensions with a circularly-symmetric trap. The GP equation describes the properties of a 
dilute trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in imaginary time over small time 
steps. The method yields solutions of stationary problems. 

New version program summary (v) 

Program title: imagtime3d 

Title of electronic files: imagtime3d.c, imagtime3d.h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 1.2 GByte 

Programming language used: C 

Typical running time: 1.5 hours (CPU1), 1 hour (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with an anisotropic trap. The GP equation describes the properties of a dilute 
trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in imaginary time over small time 
steps. The method yields solutions of stationary problems. 

New version program summary (vi) 

Program title: imagtime3d-th 

Title of electronic files: imagtime3d-th.c, imagtime3d-th.h 
Computer: Any modern computer with C language compiler installed 
Maximum RAM memory: 1.2 GByte 
Programming language used: C/OpenMP 

Typical running time: 25 minutes (CPU1), 10 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with an anisotropic trap. The GP equation describes the properties of a dilute 
trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in imaginary time over small time 
steps. The method yields solutions of stationary problems. 

New version program summary (vii) 

Program title: imagtimeaxial 

Title of electronic files: imagtimeaxial. c, imagtimeaxial. h 
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Computer: Any modern computer with C language compiler installed 
Maximum RAM memory: 32 MByte 
Programming language used: C 

Typical running time: 30 minutes (CPU1), 20 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with an axially-symmetric trap. The GP equation describes the properties of a 
dilute trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in imaginary time over small time 
steps. The method yields solutions of stationary problems. 

New version program summary (viii) 

Program title: imagtimeaxial-th 

Title of electronic files: imagtimeaxial-th. c, imagtimeaxial-th. h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 32 MByte 

Programming language used: C/OpenMP 

Typical running time: 10 minutes (CPU1), 5 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with an axially-symmetric trap. The GP equation describes the properties of a 
dilute trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in imaginary time over small time 
steps. The method yields solutions of stationary problems. 

New version program summary (ix) 

Program title: imagtimesph 

Title of electronic files: imagtimesph. c, imagtimesph. h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 2.5 MByte 

Programming language used: C 

Typical running time: 2 minutes (CPU1), 1 minute (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with a spherically-symmetric trap. The GP equation describes the properties of 
a dilute trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in imaginary time over small time 
steps. The method yields solutions of stationary problems. 

New version program summary (x) 

Program title: realtimeld 

Title of electronic files: realtimeld. c, realtimeld. h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 4 MByte 

Programming language used: C 

Typical running time: 15 minutes (CPU1), 10 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in one space dimension with a harmonic trap. The GP equation describes the properties of a dilute trapped 
Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in real time over small time steps. 
The method yields solutions of stationary and non-stationary problems. 
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New version program summary (xi) 

Program title: realtime2d 

Title of electronic files: realtime2d.c, realtime2d.h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 8 MByte 

Programming language used: C 

Typical running time: 15 minutes (CPU1), 10 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in two space dimensions with an anisotropic trap. The GP equation describes the properties of a dilute 
trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in real time over small time steps. 
The method yields solutions of stationary and non-stationary problems. 

New version program summary (xii) 

Program title: realtime2d-th 

Title of electronic files: realtime2d-th.c, realtime2d-th.h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 8 MByte 

Programming language used: C/OpenMP 

Typical running time: 5 minutes (CPU1), 2 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in two space dimensions with an anisotropic trap. The GP equation describes the properties of a dilute 
trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in real time over small time steps. 
The method yields solutions of stationary and non-stationary problems. 

New version program summary (xiii) 

Program title: realtimecir 

Title of electronic files: realtimecir. c, realtimecir. h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 3 MByte 

Programming language used: C 

Typical running time: 5 minutes (CPU1), 5 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in two space dimensions with a circularly-symmetric trap. The GP equation describes the properties of a 
dilute trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in real time over small time steps. 
The method yields solutions of stationary and non-stationary problems. 

New version program summary (xiv) 

Program title: realtime3d 

Title of electronic files: realtime3d.c, realtime3d.h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 700 MByte 

Programming language used: C 

Typical running time: 15 hours (CPU1), 12 hours (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with an anisotropic trap. The GP equation describes the properties of a dilute 
trapped Bose-Einstein condensate. 
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Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in real time over small time steps. 
The method yields solutions of stationary and non-stationary problems. 

New version program summary (xv) 

Program title: realtime3d-th 

Title of electronic files: realtime3d-th.c, realtime3d-th.h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 700 MByte 

Programming language used: C/OpenMP 

Typical running time: 4 hours (CPU1), 1.8 hours (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with an anisotropic trap. The GP equation describes the properties of a dilute 
trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in real time over small time steps. 
The method yields solutions of stationary and non-stationary problems. 

New version program summary (xvi) 

Program title: realtimeaxial 

Title of electronic files: realtimeaxial. c, realtimeaxial. h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 4 MByte 

Programming language used: C 

Typical running time: 10 minutes (CPU1), 5 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with an axially-symmetric trap. The GP equation describes the properties of a 
dilute trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in real time over small time steps. 
The method yields solutions of stationary and non-stationary problems. 

New version program summary (xvii) 

Program title: realtimeaxial-th 

Title of electronic files: realtimeaxial-th. c, realtimeaxial-th. h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 4 MByte 

Programming language used: C/OpenMP 

Typical running time: 5 minutes (CPU1), 1 minute (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with an axially-symmetric trap. The GP equation describes the properties of a 
dilute trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in real time over small time steps. 
The method yields solutions of stationary and non-stationary problems. 

New version program summary (xviii) 

Program title: realtimesph 

Title of electronic files: realtimesph. c, realtimesph. h 

Computer: Any modern computer with C language compiler installed 

Maximum RAM memory: 2.5 MByte 

Programming language used: C 
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Typical running time: 5 minutes (CPU1), 5 minutes (CPU2) 

Nature of physical problem: This program is designed to solve the time-dependent GP nonlinear partial differential 
equation in three space dimensions with a spherically-symmetric trap. The GP equation describes the properties of 
a dilute trapped Bose-Einstein condensate. 

Method of solution: The time-dependent GP equation is solved by the split-step Crank-Nicolson method by discretiz- 
ing in space and time. The discretized equation is then solved by propagation in real time over small time steps. 
The method yields solutions of stationary and non-stationary problems. 
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